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Abstract. We discuss the role of the so-called Darwin-Foldy term in the evaluation of the proton and 
deuteron charge radii from atomic hydrogen spectroscopy and nuclear scattering data. The question of 
whether this term should be included or excluded from the nuclear radius has been controversially dis- 
cussed in the literature. We attempt to clarify which literature values correspond to which conventions. A 
detailed discussion of the conventions appears useful because a recent experiment [R. Pohl et ai, Nature 
466, 213 (2010)] has indicated that there is a discrepancy between the proton charge radii inferred from 
ordinary ("electronic") atomic hydrogen and muonic hydrogen. We also investigate the role of quantum 
electrodynamic radiative corrections in the determination of nuclear radii from scattering data, and propose 
a definition of the nuclear self energy which is compatible with the subtraction of the radiative corrections 
in scattering experiments. 
PACS: 14.20.Dh, 13.40.-f, 21.00.00, 31.30.jf 
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1 Introduction 

Quite surprisingly, the recent muonic hydrogen Lamb shift 
experiment Q] at PSI has led to a value of the proton 
charge radius which is in disagreement with both the 
2006 CODATA value of the mean-square proton charge 
radius [5] , as well as in disagreement with the mean-square 
charge radius derived from electron-proton scattering ex- 
periments [3HS]. This disagreement raises a number of 
questions, two of which are the following. 

(i) Is it possible that different conventions have been 
used in order to infer the mean-square charge radius of 
the proton in atomic and nuclear physics? In particu- 
lar, in Ref. |BJ, the authors advocate to add the so-called 
Darwin-Foldy correction to the proton charge radius. Yet, 
its inclusion either into the nuclear radius [BJ or into the 
electron binding energy [7] has been the subject of discus- 
sions and is known to depend on the spin of the nucleus [7] . 
The question is whether the inclusion or exclusion of the 
Darwin-Foldy correction has been implemented consis- 
tently in all determinations of the proton charge radius 
in atomic and nuclear physics. 

(ii) The quantum electrodynamic (QED) radiative cor- 
rections to electron-proton scattering have been discussed 
in a number of papers, notably Refs. [8HT2"]. The QED 
corrections are subtracted before the form factors of the 
proton arc deduced from experiment. The notion is that 
the electric and magnetic Sachs form factors of the proton 
(Ge and Gm) should be defined so that they correspond 
to the internal structure of the proton. The same applies 
to the mean-square proton charge radius, which is pro- 
portional to the slope of the Ge form factor at zero mo- 
mentum transfer. The corresponding correction in atomic 



physics is the so-called nuclear self-energy [T3] . The ques- 
tion is whether the subtraction of the radiative correc- 
tions in scattering experiments are compatible with the 
common definition of the nuclear self-energy used in the 
atomic physics literature. 

Here, we attempt to answer both of these questions, 
and we also investigate the role of QED radiative correc- 
tions in the determination of the nuclear radius. Strictly 
speaking, the slope of the Dirac form factor F\ and of 
the Sachs form factor of the proton is known to be 
infrared divergent for any spin of the nucleus |14) . unless 
quantum electrodynamic (QED) radiative corrections are 
subtracted. From the atomic physics point of view, this 
divergence is manifest in a logarithmic term in the nu- 
clear self-energy which contributes to the atomic binding 
energy. Radiative corrections and infrared bremsstrahlung 
effects (which depend on the acceptance range of the de- 
tectors) are subtracted before evaluating the slope of the 
form factors. As a cultural matter, this aspect is not men- 
tioned in the pertinent literature [3]-[5]. 

From scattering experiments, we have a rather old 
value from Ref. |15j for the root-mean-squarc proton 
charge radius r p = yl (r 2 ) p , which reads r p = 0.88(3) fm. 
It was obtained using a dipole fit to the form factor. In 
Ref. [3J, this value has been confirmed and improved to 
yield r p = 0.880(15) fm. Then, according to Refs. @JE], 
the proton radius inferred from the world scattering data 
reads 

r p = 0.895(18) fm, (1) 
in good agreement with the 2008 CODATA value of 

r p = 0.8768(69) fm. 



(2) 



2 
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The latter value is mainly inferred from the analysis 
of spectroscopic data from atomic hydrogen and deu- 
terium spectroscopy [5J[TS] . The most recent and accurate 
measurement of the proton radius from electron scatter- 
ing [T7] , yields a value of 



0.879(8) fm, 



(3) 



when the statistical and systematic uncertainties given in 
Rcf. [17] are added quadratically, in excellent agreement 
with the CODATA value © inferred mainly from atomic 
spectroscopy and the value (P) inferred from the world av- 
erage of scattering data. However, there is a large discrep- 
ancy with the recent value from the PSI muonic hydrogen 
experiment, which reads 



r p = 0.84184(67) fm. 



(4) 



Because of this discrepancy, a study of the conventions 
used in the determination of the proton radius from scat- 
tering data and spectroscopy is indicated. 

We proceed as follows. First, the role of the Darwin- 
Foldy correction in atomic and nuclear physics determi- 
nations of the proton charge radius is analyzed (Sec. [5]). 
Radiative corrections to the proton line are discussed in 
Sec. [31 Conclusions are reserved for Sec. 2] SI (MKSA) 
units are used throughout the article unless stated other- 
wise. 



2 Darwin-Foldy Correction 

In Ref. [7J, it has been shown that the zitterbcwegung 
term of the nucleus is absent in the atomic binding en- 
ergy for spin-0 and spin-1 nuclei such as the deuteron. We 
here use the conventions of Ref. [7J, acknowledging that 
others exist [18] ■ From the nuclear physics side [5], it has 
been recommended to include the so-called Darwin-Foldy 
correction in the value of the proton radius. There is a 
connection between these two statements which will be 
explored in the following. 

The zitterbcwegung term of the nucleus is part of the 
so-called Barker-Glover corrections to atomic energy lev- 
els [15] . The Barker-Glover corrections follow from the 
two-body Breit Hamiltonian (Chap. 83 of Ref. [3U]) and 
are listed in the last term on the right-hand side of Eq. (10) 
of Ref. 0. They read 



[Zafml c 2 
2ri i m 2 N 



1 



j + 1/2 i+l/2 



(1 - 6 eo ) , (5) 



where Z is the nuclear charge number, a is the fine- 
structure constant, m r is the reduced mass of the system, 
mjv is the nuclear mass, c is the speed of light, n the 
main quantum number, and j and i are the total and the 
orbital angular momentum quantum numbers of the ref- 
erence state. Indeed, the last term in this expression is the 



Darwin-Foldy term, 



{Zafmlc 2 
2n 3 m 2 N 



E 



DF 



1 



1 



i + 1/2 £+1/2 



[Zafm lc 2 
2n 3 r 



(0) 



'N 



which is due to the zitterbewegung term of the nucleus. 
Alternatively, it can be written as 



-Edf = - I — 



m r \ (Z'a) 4 m e c 2 



i 3 Xi 



m 2 



Arrive 2 



(7) 



Here, a is the fine-structure constant, and Z is the nuclear 
charge number, m r is the reduced mass of the system, m e 
is the electron mass, c is the speed of light, and Ac is the 
Compton wavelength of the electron divided by a factor 
2tt. The Darwin-Foldy term is nonvanishing only for S 
states (£ = 0). 

The main nuclear-size shift of atomic energy levels is 
given by the well-known expression 



2 / m r 

-C'NS = o 

3 \m e 



(Za)^m e (? 
n 3 X 2 n 



{r 2 ) N S m , 



(8) 



where (r 2 ) n is the mean-square charge radius of the nu- 
cleus. The expressions ([7]) and ([5]) have the same structure. 
The nuclear size effect and the Darwin-Foldy correction 
have the same effect on the spectrum if we alternatively 
add the Darwin-Foldy correction to the atomic energy 
levels or if we add the term in curly brackets in Eq. ([7|) 
to the mean square nuclear radius. Let us refer to the nu- 
clear radius defined without the Darwin-Foldy term as the 
nuclear radius in "atomic physics" (ATP) conventions. In- 
deed, this definition has implicitly been proposed in the 
paper |21j . where the nuclear radius difference of proton 
and deuteron was experimentally measured and theoreti- 
cally evaluated. The following two replacements are found 
to be equivalent, 



-Ens — > -Ens 
where 



Ebf & (r 2 ) N 



ATP 



DF 



• 2V DF 3ft 

4m^c 2 



(9) 
(10) 



Here, N denotes the nucleus (N = p for the proton), and 
rriN is the nuclear mass. According to Ref. [6], in atomic 
physics conventions, the proton charge radius is propor- 
tional to the slope of a subtracted Sachs form factor Ge, 



„2\ATP 



(r 



2\p 
E 



dG E {q 2 ) 



dq 2 



(11) 



<j 2 =0 



where q 2 = (q ) 2 — Q 2 is the momentum transfer. In or- 
der to be consistent, it is important to stress that a sub- 
tracted form factor Ge has to be used in the evaluation of 
the slope, because formally, the Sachs form factor Ge has 
an infinite slope at zero momentum transfer, due to an 
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infrared divergence at the one-loop level, which is caused 
by vertex corrections involving virtual interactions with 
very soft virtual photons. This is illustrated in the dis- 
cussion following Eq. (|15|) below. The infrared divergence 
thus is of quantum electrodynamic origin and not a con- 
sequence of the internal structure of the proton, and it 
only enters at the one- loop level (order a). In the litera- 
ture, the subtractions are sometimes carried out tacitly, 
and it is understood that the form factor G E employed 
for the proton is due entirely to its internal structure, and 
all QED effects have been subtracted. 

Let us set this problem aside for the moment and con- 
tinue to study the role of the Darwin-Foldy correction in 
the definition of the nuclear charge radius. An alterna- 
tive convention for the proton charge radius is being dis- 
cussed in Ref. [BJ • We refer to this convention as the Friar- 
Martorell-Sprung (FMS) convention as it has been advo- 
cated in Ref. [BJ. In FMS conventions [see also Eq. (O], the 
mean square proton charge radius includes the Darwin- 
Foldy correction, 



/ 2 >FMS / 2 \ ATP / 2\ DF 

( r ) P = ( r ) P + ( r ) P 



(12) 



According to Ref. [BJ, the alternative mean square charge 
radius can be written as the slope of a modified Sachs form 
factor 



G E (q 2 ) 



G E (q 2 ) 



^l-q 2 /Am 2 p 



6h 



2 dG E {q 2 ) 



dq 2 



(13a) 



(13b) 



The authors of Ref. [BJ show that this representation is 
better adapted to the relativistically invariant representa- 
tion of the Rosenbluth [22] formula. 

If one uses the FMS convention (TT2)) for the proton 
charge radius, then, even for a pointlike nucleus, the nu- 
clear size correction is nonvanishing for the atomic bind- 
ing energy. Therefore, the ATP definition of the charge 
radii has been favored in Ref. |23| . who have argued in 
Sec. 7.1.1 of the cited literature reference that in Ref. [BJ, 
"it is suggested to include the Darwin-Foldy contribution 
in the definition of the nuclear charge radius. While one 
can use any consistent definition of the nuclear charge ra- 
dius, this particular choice seems to us to be unattrac- 
tive since in this case even a truly pointlikc particle in the 
sense of quantum field theory (say an electron) would have 
a finite charge radius even in zero-order approximation." 
Here, the authors refer to the zcroth-order approximation 
as the one without any radiative corrections. In order to 
ease our mind, we may note that the Darwin-Foldy cor- 
rection vanishes for an infinitely heavy, point nucleus, i.e., 
in the limit of — > oo. In this limit, the ATP and 
FMS conventions for the charge radius definition are in 
agreement. Due to the uncertainty principle, one cannot 
locate a particle and therefore, its charge, to better than 
its Compton wavelength. The Darwin-Foldy correction is 
of this magnitude. 



There is thus an immediate need to clarify which con- 
ventions have actually been used in the determinations of 
the proton charge radius in Refs. [5HS| and in the CO- 
DATA adjustment in Ref. [2J. We observe: 

(i) According to the first (unnumbered) equation given 
on p. 410 of Ref. [BJ, one may infer that in Refs. 3 5 , the 
form factor G E (or merely an infrared safe, subtracted 
variant of it) and not G E is being used for the determi- 
nation of the rms radius given in Eq. ([1}; this observa- 
tion can be confirmed [21]. This means that the values 
given in Refs. J3HS] and in particular in Eq. (fTJ) are in 
agreement with the atomic physics conventions mentioned 
above. The same applies to the values given in Eq. ([3]). 

(ii) In the latest CODATA adjustment given in Ref. [2], 
the conventions are consistent with those used here and 
in Ref. [21], although this is somewhat less obvious. First 
of all, we reemphasize that, as evident from Eq. (10) of 
Ref. [2J, the atomic binding energy is defined to include 
the Dirac-Foldy term in the CODATA adjustment. For 
the deuteron, there is no Darwin-Foldy correction in the 
atomic binding energy (see [7]). In Ref. [2J, the authors still 
use the Darwin-Foldy correction in the atomic part of the 
energy (even for deuterium) and later add the Darwin- 
Foldy correction back on to the deuterium radius. This is 
consistently done in all CODATA adjustments since the 
1998 adjustment (see Ref. [2S]), and a pertinent discussion 
can be found in Appendix A8 near Eq. ( A56) of Ref. [2BJ . 
A superficial reading of Appendix A8 of Ref. [25] might 
otherwise suggest that the Darwin-Foldy correction has 
been inadvertently added to the deuteron (rather than 
proton!) nuclear charge radius in Ref. [2BJ; but this addi- 
tion is compensated by the inclusion of the Darwin-Foldy 
term even for bound-state deuterium energy levels, where 
according to Ref. [7J, this term should have been excluded. 
Therefore, both proton and deuteron charge radii given in 
Ref. [2] correspond to ATP conventions [21] and exclude 
the Darwin-Foldy term from the nuclear radii for both 
proton and deuteron. 



3 Radiative Corrections 

We now turn our attention to the role of radiative correc- 
tions to the proton line in both nuclear physics scatter- 
ing experiments as well as in the determination of atomic 
binding energy levels. We recall once more that in atomic 
physics (ATP) conventions, the proton charge radius is 
defined as the slope of a subtracted electric G E Sachs 
form factor, (r 2 )£ TP = 6h 2 dG E /dq 2 \ q 2 =0 . The electric 
and magnetic G E and Gm form factors for the proton 
(a spin-i particle) are related to the Dirac F\ and Pauli 
Fi form factors by the relations 



^ 2 )=i^ 2 ) + i ^W). 



(14a) 



G M (q 2 ) = F 1 (q 2 ) + F 2 {q 2 ), F 2 (0) = x p , (14b) 

where x p = (g p - 2)/2 = 1.792 847 356(23) gives the 
anomalous magnetic moment of the proton (see Ref. [2]). 
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Fig. 1. Radiative vertex correction to electron-proton scattering (heavy line) and associated 
brcmsstrahlung diagrams. 



For a point particle like the electron, we have x e = 
a/(2n) -C 1, but for a particle like the proton, the bulk 
of the contribution to x v is from the internal structure, 
whereas a tiny correction also is due to its electromag- 
netic nature. 

Let us therefore make the following separation, 



G E (q 2 ) = G E (q 2 ) + Gf B (q 2 ), 
^i(9 2 )=^i(9 2 ) + ^ QED (9 2 ), 

^ 2 ) = iW) + ^ QE V)- 



(15a) 
(15b) 
(15c) 
(15d) 



Here, the ovcrlined quantities represent the contributions 
to the proton form factor due to its internal structure, 
whereas the quantities with the superscript QED repre- 
sent the contributions due to QED point-particle theory. 

The one-loop slope of G§ ED is infrared divergent solely 
due to QED vertex corrections (sec Fig. [1}, but the slope 
of Ge is infrared safe. Indeed, the slope of the Dirac iq 
form factor of any charged spin-^ particle is infrared di- 
vergent, and the coefficient of the logarithmic infrared di- 
vergence of F\ even is independent of the nuclear spin |14j . 
The infrared problem persists both in the atomic physics 
determination of the nuclear charge radius as well as in 
the determination from scattering data. From Chap. 7 of 
Ref . [13] , we know that the pure QED radiative contribu- 
tion from the first diagram in Fig. [T] corresponds to the 
following replacement of the Dirac current 7 M of the pro- 
ton, 



Y + < 7 



^ QED V) - 1 



2m r 



-F, 



QED 



(16) 

where F-P EE> and F 2 ^ ED are the QED expressions for the 
Dirac and Pauli form factors of a spin-i point particle 
with the proton mass, as given in Eqs. (7-60) and (7-58) 
of Ref. [13], respectively. In particular, we haveF 2 QED (0) = 
a/(2ir). From Eqs. (3.36) and (4.14) of Ref. [9], and from 
Eq. ( A76) of Ref. jTD] , we may conclude that the radiative 
QED term in curly brackets in Eq. (fTo]) is indeed sub- 
tracted when the radiative corrections to electron-proton 
scattering arc eliminated from experimental scattering 
data. 

We now have to turn our attention to atomic physics 
and identify the nuclear self-energy as the correction to 
atomic energy levels corresponding to the terms in curly 



brackets in Eq. (fTl)]) , which arc subtracted from scattering 
data. To this end, we first make a slight detour and observe 
that the proton radius definition according to 



„2\ATP 



= 67r 



6h 



dG E (g 2 ) 
dq 2 

2 dF^q 2 ) 



dq 2 



q 2 =0 



q 2 =0 



m 2 



2 (m p c) 



2 ' 



(17) 



entails a term originating from the "internal" contribution 
Tip to the anomalous magnetic moment of the proton. We 
note that vc v « Tc p because x^ ED = a/ (2n) <C x v is small. 
The anomalous magnetic moment of the electron shifts S 
state energy levels [see Eqs. (42) and (44) of Ref. [53]] by 



a \ (Za) 4 



'•"r \ 2 r 



(Za) 4 



in i 
in, 



m e c 2 S, 



m ■ 



(18) 



The corresponding proton line contribution reads 



2 r 

m r c Og, 



2 (m r \ 3 {Za) A m e c 2 /3 ft 2 _ 
~3\m^J ti3 X 2 C \2 (m p c) 2 Kp 



(19) 



The slope of the Fi form factor of the proton leads to an 
energy correction, 



E E , — 



2 (m r \ (Za) 4 m e c 2 



n*\ 2 c 



2 dF^q 2 ) 



dq 2 



q 2 =Q 




The sum of E x and E El is 
3 



2 / TO r ^ (Za) 4 m e c 2 

3 \ m» 




(21) 
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which corresponds to the well-known expression @ and 
clarifies that indeed, the anomalous magnetic term of the 
proton forms part of the nuclear size correction in the 
ATP conventions, and that indeed, this convention is in 
agreement with the mean-square-radius of the proton be- 
ing defined as the slope of the Sachs Ge form factor. The 
expectation value of the Dirac S for S states is 



ir(Za)S 3 (r] 



(Za) 



nS 



m r c 



(22) 



In order to describe the two-body interaction includ- 
ing the nuclear-size effect, we have to consult the two- 
body Breit Hamiltonian, and temporarily switch to nat- 
ural units with ft = c = e$ = 1. The Schrodinger Hamil- 
tonian of the two-body system consisting of an orbiting 
spin-i particle of mass m e (the electron) and a spin-^ 
nucleus of mass m^v and nuclear charge number Z is 



p 2 Za 

H S = ■z , 

lm r r 



(23) 



where m r is the reduced mass of the system. The two- 
body Brcit-Pauli Hamiltonian with anomalous magnetic 
moment corrections for electron and proton reads 



H = 



4 4 
p p 



8ml 8m 3 
2Za ( 3 3 



/ 3 3 , 9n atp\ -o. . 



Za (& 3 r l ri 



2m e m p r r I r ' r 3 



P — 



Za 



Hyperfine structure effects are by definition excluded from 
the Lamb shift [26] . 

We now have to carefully define the nuclear self-energy 
for S states so that the QED contributions to the proton 
factors, which were excluded from the form factors used in 
Eq. (|T9|) . (|20|) . and (|2"Tj) for the evaluation of the nuclear- 
size correction, are included instead into the nuclear self- 
energy. Furthermore, we notice that the "fs" terms, which 
are proportional to the spin-orbit coupling L ■ cr e of the 
electron, arc defined with the phcnomcnologically inserted 
anomalous magnetic moment >c e of the electron. We also 
have to define the nuclear self-energy so that the terms 
due to the >r^ ED = F 2 ^ ED (0) part of proton anomalous 
moment, which are already included in the phenomeno- 
logically inserted full anomalous magnetic moment x p in 
the L ■ cr p terms in Eq. (f2"4"]) , are not double counted in 
the nuclear self-energy. 

The F\ form factor of the electron induces the following 
correction to the Lamb shift, 



E R = 



a{Za) 4 



J 4 lnf^ 
V 2e 



11 

18 



1 



(25) 

where e is a noncovariant, infrared cutoff parameter [271 
|2"8"] . The matching with the covariant photon mass has 
given rise to some discussion |29j in the early days of quan- 
tum electrodynamics. The terms 11/18 and 1/2 are from 
the nonlogarithmic term of the slope of the F\ = F® EU 
form factor of the electron, and the term 1/2 is due to 
From the proton line, we have the following, 
corresponding contribution due to the QED parts F^ EU 
and F® EU of the proton form factors, 



Za Za 

+ 0- + _3 £ ' ff « + (! + 2X P)Z^3 £ • a P 



2m e m p r 3 



Arn^r^ 



ts hfs 
+ (1 + X p )- • <T p H 

im^.rripr'' Amempr^ 



his 



his 



, 8tt 3 <r e • r cr p ■ r (T e cr p 
X {—<T e <Tp 6 3 (r)+3 ^ ^ 



(24) 



Ill's 



Here, we keep the nuclear charge number Z as a variable 
so that the expression below can be readily generalized 
to a spin-i nucleus with Z / 1. Terms labeled with "fs" 
are relevant for the fine structure, whereas terms labeled 
by "hfs" correspond to the hyperfine structure. All terms 
in Eq. (IM1) proportional to the electron anomaly x e con- 
tribute to the fine structure and thus, to the Lamb shift of 
these states. Except for the term incorporated into the nu- 
clear mean-square charge radius (f 2 ) f) , the terms pro- 
portional to the (complete) proton anomalous magnetic 
moment >c p listed in Eq. (|24[) influence only the exchange 
of a magnetic photon of the nucleus and the orbiting par- 
ticle, i.e., the hyperfine structure (not the Lamb shift). 



E 



NSE.l = 



Z{Zaf 



m r c 



4 ln pA 10 
3 V 2e J 9 



(26) 

for S states. The infrared divergence in (|2"6"|) . which for 
scattering experiments is compensated by bremsstrahlung 
diagrams, is cut off in atomic physics at the binding energy 
scale {Za) 2 m r . The matching low-energy contribution |30j 
contains the Bcthc logarithm hik^n^t). 



E 



NSE.2 



Z(Zaf 



2e 



m r c 



(27) 



The sum of ((26]) and ((27]) is free from the scale separation 
parameter e and reads 



Z(Zaf f m r \ 2 

£nse = ^— I rn r c 

Trn a \ nip / 



(28) 



hi 



m r (Za)' 



10 



Seo - g lnfc (n,^) } . 



This expression generalizes the nuclear self-energy to 
states with nonvanishing angular momenta. 
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A few remarks arc in order. First, we observe that the 
anomalous magnetic moment x e , due to the spin-orbit 
coupling of the electron, leads to the following term in 
the self-energy for non-5 states, 

„, afZa) 4 ( m r \ 2 , / 1 \ . „ , 

E = "W" — m * c - o (opm^ 1 ~ Seo) ' 

nn J \m e ) \ 2k(2i + 1)J 

(29) 

where k — 2(jE—j) (j + 1/2) is the Dirac quantum number. 
The analogous term due to the proton line is given in 
Eq. (153) of Ref. [121, 

However, this term originates from the spin-orbit coupling 
of the proton and is already contained in the phenomeno- 
logically inserted x p in Eq. (j24]) . We therefore exclude this 
term from the definition of the nuclear self-energy for nan- 
s' states. 

The Bcthc logarithm term lnfco(n, t) in Eq. ([28]). for 
a general hydrogenic state, can be derived on the basis of 
nonrclativistic perturbation theory for the two-body sys- 
tem [31 j - It is a consequence of the low-energy part of 
the proton self energy, which is due to the exchange of 
low-energy photons along the proton line. The logarith- 
mic term (in e) in the low-energy part given in Eq. (|27p 

compensates the infrared divergence of the F 1 QED form 
factors at the atomic binding energy scale, in analogy 
to the cancellation of infrared divergences in the vertex 
and bremsstrahlung diagrams (for scattering). The defini- 
tion (1281) takes into account the quantum electrodynamic 
properties of the proton (as a charged spin-^ particle) 
and is compatible with the radiative corrections listed in 
Eq. (JTSJ). The nonlogarithmic term 10/9 in Eq. (|2g]). which 
constitutes an addition to expression for the nuclear self- 
energy given in Ref. [T3], corresponds to a shift of the 
proton radius by about 5r p = 0.0001 fm or alternatively 
to 0.20 kHz in frequency units for the hydrogen-deuterium 
isotope shift of the 1S-2S transition. 

Finally, one might argue with respect to our result (|28[1 
that the starting point of QED calculations connected 
with the proton is a proton with form factors that are 
due to strong interactions, and that, therefore, all QED 
integrals are cut off from above by the proton radius, i.e., 
that the proton structure significantly influences the ex- 
change of virtual photons of wavelengths shorter than its 
own size. However, the e parameter is an infrared cutoff, 
and it stems from the infrared divergence of the if ED 
form factor of the proton. Indeed, the physically appro- 
priate values for the e parameter are in the region 

{Za) 2 m e <Cf<m e <i, m p , (31) 



where A = VOJlGeV 2 = 0.842GeV is a parameter char- 
acterizing the proton size. The first inequality in Eq. (|31[) 
is due to the requirement that all bound-state poles must 
be integrated over in evaluating the low-energy part of 
the self-energy, and the latter inequality is due to the 



rclativistic nature of the electron mass scale (separation 
from the high-energy part). The coefficient multiplying 
the logarithmic infrared divergence in Eq. (f21)]) is the 
same whether we use ln[yl/(2e)] for the logarithmic term, 
as done in Eq. (154) of Ref. [23], or ln[m p /(2e)], with 
m p = 0.938 GeV, as done here. The proton structure thus 
docs not influence the leading coefficient of the logarith- 
mic divergence of the infrared behavior of the form factor. 
The question then is which portion of the proton struc- 
ture should be counted as a radiative nuclear self energy, 
and which should be counted as a proton structure cor- 
rection. Here, we advocate the view that if the radiative 
corrections to scattering are accounted for by the if ED 
and if ED contributions to the form factors of the pro- 
ton (treated as a point-like QED particle), then, also, the 
proton should be treated as a point-like QED particle in 
evaluating the nuclear self energy in atomic physics. The 
difference of the above mentioned two logarithms, 

Ks) --(£)--(£)■ (32) 

is infrared safe, part of F\ and therefore, part of the 
nuclear-size correction in our conventions. While this view 
is at variance with the result obtained in Eq. (154) 
Ref. [23], we stress that the precise definition of the nu- 
clear self energy depends on the way in which radiative 
corrections to the proton line are subtracted in atomic 
physics and in scattering experiments. Our view thus holds 
relative to the conventions used for subtracting radiative 
corrections to the proton line in scattering experiments. 

In the evaluation of higher-order nuclear structure cor- 
rections to the bound-state spectrum, wc only have to 
remember that the Ge inferred from scattering experi- 
ments is not the full Ge of the proton, but with the con- 
tributions from if ED and if ED subtracted. The QED 
contributions to the proton form factors thus have to be 
treated separately. One example can immediately be given 
as follows. It is known [32] that the third Zcmach moment 
correction [33J to bound-state energy levels is obtained in 
second order perturbation theory involving the correction 
to the Coulomb potential due to the finite size of the nu- 
cleus, or in other words, in the second order with respect 
to the slope of the proton form factor. The contribution 
from the radiative correction due to the QED contribu- 
tion F-P ED (g 2 ) to the Dirac form factor of the proton is 
thus neglected in the derivation of the Zcmach correction. 
If desired, this effect can be added back on. It is easy to 
write down the leading logarithm 

bW( „ s)= » „j 

27 tt z n 6 \Tn>pJ 

for the energy correction to S states in the second order of 
perturbation theory with respect to the proton self-energy. 
The numerical value of this correction is tiny because it is 
of high order in the electron-proton mass ratio. 
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4 Conclusions 

In this paper, we have analyzed the role of the Darwin- 
Foldy term in the determination of atomic energy levels 
and in nuclear physics scattering experiments that mea- 
sure the proton form factors. In Sec. [5J we have clarified 
that all current determinations of the proton charge ra- 
dius, both in atomic as well as in nuclear physics, are based 
on the slope of the electric Sachs Ge form factor (jlip and 
thus exclude the Darwin-Foldy term. This applies to both 
the value given in Eq. ([I} as well as the value given in 
Eq. §2$. The latter value is determined on the basis of 
a comparison of the most accurately measured transition 
frequencies in atomic hydrogen to theory [16] . 

Furthermore, we have clarified that although a su- 
perficial reading of Appendix A8 of Ref. [22] might sug- 
gest that the Darwin-Foldy correction has been inadver- 
tently added to the deuteron (rather than proton!) nuclear 
charge radius in Ref. [35] , this is actually not the case: the 
Darwin-Foldy correction is absent for deuterium [7], and 
if one includes this contribution into atomic energy levels, 
even for deuterium, as done by the authors of Refs. [21I25]. 
then one has to add the Darwin-Foldy correction back 
on to the deuteron radius as is done in Appendix A8 of 
Ref. [25] (as well as, consistently, in all CODATA adjust- 
ments since 1998). 

These observations imply that the Darwin-Foldy term 
is excluded from both values of the proton radii given in 
Eqs. (TIJ and ([5]), which in addition are in good mutual nu- 
merical agreement. In Ref. [T], the proton radius is given 
as r p = 0.84184(67) fm, in disagreement with Eqs. (TIJ 
and @, and with the Darwin-Foldy term likewise being 
excluded. A conceivable accidental incompatibility of the 
conventions used in Refs. [THE] for the proton radius there- 
fore cannot be the reason for the observed discrepancy. 

In Sec. [3] we turn our attention to radiative corrections 
to electron-proton scattering. Indeed, the infrared diver- 
gence of the Sachs Ge form factor at the one- loop level en- 
ters both the atomic physics determination of the proton 
radius as well as the nuclear physics determination from, 
e.g., scattering data. While it is clearly stated in Refs. 3 5 
that the proton charge radius is determined based on scat- 
tering data for Ge, it is useful to remark that the slope 
of the Ge form factor is determined in each case from 
subtracted scattering data: namely, the infrared divergent, 
QED radiative corrections whose exact values depend on 
the acceptance range of the detectors are subtracted when 
the slope of the form factor is being calculated. Theoretical 
input data for the radiative corrections is available from 
a number of sources, notably, Refs. [8HT2]. Equation ((28]) 
represents an attempt to define the nuclear self energy in 
full compatibility with the subtractions of radiative cor- 
rections to electron-proton scattering, notably, Eqs. (3.36) 
and (4.14) of Ref. 0, and Eq. (A76) of Ref. [10]. The re- 
sult ([2"5]) also represents a generalization of the nuclear self 
energy to non-S 1 states (i ^ 0). 

Finally, let us remember that the nuclear self energy is 
a numerically rather tiny correction, especially for systems 
with a small mass ratio of the orbiting particle versus the 
nucleus. Likewise, the consistent subtraction of the radia- 



tive corrections listed in Eq. (|16[) in a scattering experi- 
ment depends on the availability of accurate data at low 
momentum transfer, and on the accurate subtraction of 
other contributions, such as double scattering (including 
effects due to the polarizability of the proton) and double 
bremsstrahlung. It is known [SHU] that the contribution 
of the latter effects may be numerically more significant 
than that of the QED terms in curly brackets in Eq. (TIP]) . 
Currently, the uncertainty induced by the nonlogarithmic 
terms is negligible on the scale of the total uncertainty of 
the proton radius, but still, it is an important conceptual 
issue to carry out the analysis of radiative corrections to 
the proton line consistently in atomic and nuclear physics 
experiments. 
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